HEPHY-PUB 806/05 
|hep-ph/0504109| 



Precise predictions for chargino and neutralino pair 
production in e + e~ annihilation 



W. Oiler, H. Eberl, W. Majerotto 



Institut fiir Hochenergiephysik der Osterreichischen Akademie der Wissenschaften, 
A-1050 Vienna, Austria 



Abstract 

We present a high precision calculation of chargino and neutralino pair pro- 
duction at e + e~ colliders. Within the Minimal Supersymmetric Standard Model, 
the full one- loop and higher order QED corrections are included. Special care has 
been taken in the definition of the Lagrangian input parameters. Furthermore, the 
proper inclusion of QED corrections and the separation of weak and QED correc- 
tions are discussed. We show numerical results for total cross sections, as well as 
forward-backward and left-right asymmetries for the SPSla' scenario as proposed in 
the SPA project. The complete corrections are about 10% and in some cases even 
larger, in particular for XiX^j production with sizeable higgsino components. These 
corrections have to be taken into account in a high precision analysis. 



1 Introduction 



Supersymmetry predicts the existence of fermionic partners to the gauge and Higgs 
bosons. In the Supersymmetric Standard Model (MSSM), one has two charginos Xi 
and xt> which are the fermionic mass eigenstates of the supersymmetric partners of the 
W ± and the charged Higgs states Hf 2 - There are four neutralinos — xl being the 
fermionic partners of the photon, Z° boson, and the neutral Higgs boson states H® 2 - At 
tree-level, the chargino and neutralino systems depends only on the parameters M, M' , 
[A, and tan/? = Vi/v\, with t>i j2 the vacuum expectation values of the two neutral Higgs 
doublet fields. 

It is expected that at least some of these particles will be detected at LHC or Tevatron, 
most likely in the cascade decays of gluinos and squarks. The properties of charginos and 
neutralinos can be studied with high accuracy at a high energy linear e + e~ collider pQ. 
At lowest order, the above mentioned parameters can be extracted from the masses and 
production cross sections in e + e~ collisions with polarized beams [3 El- The experimental 
accuracy of measurements to be obtained at a linear collider will be so high that it is nec- 
essary to incorporate effects beyond leading order in the theoretical calculations in order 
to match the experimental precision. Beyond tree-level, the definition of the parameters 
involved are no more unique and depends on the renormalization scheme. Therefore, a 
well defined theoretical framework has been recently proposed within the so-called SPA 
(SUSY Parameter Analysis) project [I]. The proposed "SPA convention" provides a clear 
base for calculating masses, mixings, decay widths and production processes. The final 
goal is to extract the fundamental supersymmetry parameters from data. 
In this paper, we treat in detail the calculation of chargino and neutralino pair production 
at full one-loop level within the MSSM. Particular attention is paid to a proper inclusion 
of QED corrections and to a suitable separation of weak and QED corrections. For the 
initial state radiation of photons, we use the structure function formalism (THj, where also 
non-leading log terms are included and the soft-photon contributions are summed up to 
all orders in perturbation theory. We adopt the SPA convention, in which the SUSY input 
parameters are defined in the DR scheme at the scale Q=l TeV. The parameters are then 
translated to the on-shell scheme in which the calculations are performed. The result for 
the observables (cross sections, asymmetries, etc.) is therefore independent of the renor- 
malization scheme up to higher order corrections. The masses of the SUSY particles and 
the Higgs bosons are of course defined as physical pole masses. 

This paper is organized in the following way: Section ^introduces the used notation and 
conventions. In section |3] the tree-level cross section is given. Section 0] contains the 
calculation of virtual corrections and describes in detail the used renormalization scheme. 
In section |5] the inclusion of the QED corrections and the separation of the weak and 
QED part is discussed. Section El shows the numerical analysis for the SPSla' scenario, 
proposed within the SPA project. Finally, section [7| summarizes this paper. 



2 



2 Notations 



Y = 



In supersymmetric extensions of the SM, the fermionic superpartners of the vector and 
Higgs bosons, the gauginos and higgsinos, can mix and form the mass eigenstates of the 
neutralinos and the charginos. In the MSSM the mixing is defined by the two mass 
matrices, which are non-diagonal after the gauge symmetry breaking. The symmetric 
tree-level mass matrix of the neutral tp® = — i\ 3 , ^ , ^frj Weyl states 

/ M' —m z sin 6 W cos (3 mz sin 6 W sin (3 \ 

M mzCOsOwCOsfl — m z cos 9w sin (3 

—mz sin 6 W cos (3 m z cos 6 W cos (3 — /i 

\ m z sin 6 W sin (3 —m z cos 9w sin (3 —fj, J 

(1) 

can be diagonalized by the unitary matrix iV 

diag(m^o, m^o, m^o, m^o) = N*Y N' . (2) 
The mass eigenstates in the 4-component Majorana spinor notation Xi are defined by 

x?=^?, x*°=(sy (3) 



X 



The diagonalization of the chargino tree-level mass matrix 

x _f M V2m w cosf3\ 
m w sm(3 /j, J 

can be performed by two unitary matrices 

diag(m -± , m^±) = U* X VK (5) 

Thus, the chargino mass eigenstates in Dirac form xt is obtained from the Weyl states 
ip + = (—i\ + ,vp} l2 ), = (— ip 2 Hl ) by the relations 

Xt = V^f , XT = U^j , xt = ( f f ) ■ (6) 

The spectrum of the neutralinos and charginos is specified by the soft SU(2) and U(l) 
gaugino mass parameters M and M', the Higgs/higgsino parameter \i and tan/3. 
We assume no flavour mixing in the sfermion sector and further neglect the selectron 
L-R mixing, which is suppressed by a factor m e in the off-diagonal selectron mass matrix 
elements. Thus, we get for the selectron masses and the sneutrino mass in terms of 
electroweak parameters and SUSY breaking masses at tree-level 



m i L 

,2 



Mf i - m| cos 2(3 (\ - sin 2 9 W ) (7) 
m\ R = M| - m| cos 2(3 sin 2 9 W (8) 
ml = Ml + lml cos 2(3 (9) 
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3 Tree-level 



The tree-level pair production processes for charginos, Fig. [TJ and neutralinos, Fig. El 

e~(pi) e + (p 2 ) - xr(A:i) (hj = 1, 2) , 

e~ (Pi) - X- (*i) X° (to) (i, J = 1, 2, 3, 4) , 

have been already extensively discussed in the literature [HJ E] . After a Fierz transfor- 




Figure 1: Tree- level chargino production 




Figure 2: Tree-level neutralino production 



mation of the sfermion exchange channels, where the electron mass is neglected in the 
Yukawa couplings, the tree-level matrix-elements can be written in terms of four indepen- 
dent helicity amplitudes 

M% = i e ^Q^[v{p 2 ) 1 ^P a u{ Pl )}[u{k l )^P p v{k 2 )} , {a,(3}e{L,R}. (10) 
The bilinear charges Q a0 for charginos and Q% for neutralinos are 
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Qll — $ij + Dz Cl Uij , Qlr — + Dz Cl Vij + Do Vij , 

Qrl = $ij + Dz CnUij , Qrr = $ij + Dz Cr Vij , 



Q° LL = D z C L Uij - D U>L dj , Q° LR = -D z C L A/>- + D tjL £*■ , 

Qrl = Dz C R A/" tJ + D t>R K l3 , Q% = -D z C R N*j - D U>R 11% . 

The first index denotes the chirality of the e current, the second one of the Xi,j current. 
We have introduced the projection operators Pl/r = 1(1 "FTs)? the kinematical variables 

s = (pi+p 2 ) 2 , t=(pi-A;i) 2 and u = (p x - k 2 ) 2 , (11) 

the normalized propagators 

s s s s 

D z = n" ' D>u = , D t L i R = g — > D u L ir = 2 — j (12) 

s — m? t — m, } 2 t — mi u — m~ 

V e L/R e L/R 

and the coupling matrices 

= (S\V — g)/ {Syy Cyy) , Cr = 1/ Cyy , 

% = (sir Sij - UaU*, - \ U l2 U* 2 ) , V l3 = Uij{U -> V) , 

% = Ki^-i/(2 s^) , A/ij = (iV^ - , 
£*j = (AW + N a s w )(N* 2 c w + N* lSw )/(4 s 2 w c 2 w ) , Uij = NaNJJc^ . (13) 

Summing up the final state helicities Ai^ = Y^a=±M^p , the integrated tree-level cross 
section for polarized beams reads (for m e — > 0) 

J " a=± 

with £ ± the degrees of polarization of the e ± beams, and X(x, y, z) :— (x — y — z) 2 — <iyz. 



4 One-loop corrections 

For a high precision analysis of the neutralino and chargino sector, the inclusion of higher 
order corrections is mandatory. An appropriate regularization scheme to preserve super- 
symmetry (at least at one- loop level) is dimensional reduction (DR). In the following 
elaborate calculations, a large number of Feynman diagrams are involved. This makes 
it necessary to use an appropriate computer algebra tool. The FeynArts 3.2 and Form- 
Calc 4 |3] packages are adopted for the generation of diagrams and amplitudes. This is 
performed in the £ = 1 Feynman-t'Hooft gauge. We have further integrated our renor- 
malization scheme discussed in the next section into these programs. The computation 
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of the one- loop integrals is based on the packages LoopTools and FF [H] . 
The virtual corrections for polarized beams can be written in the form 

/ da^ = ^^^r^ / dO]P 1(1 + a£_)(l - 2Re{(M^M1} . (15) 

The one-loop matrix element Aif consists of all possible vertex corrections (Fig. Hj), self- 
energy (Fig. EI) and box (Fig. EJ) diagrams and the corresponding counter terms. We again 
neglect the electron mass m e wherever possible. 




Figure 3: Generic virtual diagrams. The virtual corrections are structured into vertex, 
propagator and box contributions. 



4.1 Renormalization scheme 
4.1.1 Gauge sector 

Since the minimal supersymmetric extension of the SM does not introduce further gauge 
couplings, the renormalization used for the SM parameters can be taken from jjj. The 
quark mixing matrix is assumed to be diagonal. The bare parameters are split into 
renormalized parameters and their counter terms. 

m 2 w ^ + 5m 2 w , m% ->• m| + 5m z , (16) 
f L ^(l + i5Z L )f L , f R ^(i + l5z R )f R , (17) 



Z \ ( 1 + UZ ZZ UZ ZA \(Z 



2 Zj Zj 2 

A) ^ \ \5Z AZ 1 + \5Z AA j \ A 



Since we assume no flavour mixing, no off-diagonal wave function counter terms for the SM 
fermions have to be introduced. The renormalization conditions that the on-shell masses 
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Figure 4: Generic diagrams for the different vertex corrections. Seven different classes of 
particles are introduced: / for all SM fermions, f\ for e and v e only, / for all sfermions, 
fx for cl.r and v e only, x f° r neutralinos and charginos, V for vector bosons, and H for 
Higgs and Goldstone bosons. 



are the real parts of the poles of the corresponding propagators and that the renormalized 
fields are properly normalized, yield the following results for the counter terms. 

5m 2 z = Renf z (mi), 5m 2 w = ReU^ w {m 2 w ) , (19) 

SZ VV = -Retl% v (m 2 v ) , V = A,Z,W, (20) 

2fenf(m|) 2 Re 11^(0) 

dZ AZ = = , dZ Z A = n ' I 21 ) 

m z m z 

with the transverse self-energies Hr and ti{m 2 ) = [^^[(fc 2 )] fc2 _ m2 . Please note that Re 
only takes the real part of the loop integrals and does not affect the possibly complex cou- 
plings. The weak mixing angle is fixed by the usual on-shell condition = 1 — rriy^/m^. 
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Figure 5: Generic self-energy diagrams for the s-channel vector and the t/u channel 
slepton propagators. FP denotes the class of Faddeev-Popov ghosts. All other particle 
classes are the same as in Fig. HJ 




Figure 6: Generic box diagrams. The notation is taken over from Fig. 0] 



S 



For the SM fermions we find the wave function counter terms 



5Z l 



Re 



n L (mJ) - m)(tl L (m)) + tl R (m 2 f )) + ^-(U SL (m 2 f ) - Tl SR (m 2 f )) 



f 



-m f (U bL (mj) + U bH (mj)) 
5Z R = 5Z L {L^R), 
where we used the decomposition 

U(k 2 ) = ]kP L Tl L (k 2 ) + ^tP R U R (k 2 ) + P L U SL (k 2 ) + P R U SR (k 2 
For the eiectric charge defined in the Thomson limit a = e 2 / (4tt 

% 



(22) 
(23) 

(24) 



- = ~SZ AA - -^5Z 
e 2 2c w 



ZA 



1/137.036, one obtains 
(25) 



However, this definition leads to large radiative corrections for processes at the GeV or 
TeV scale. Furthermore, the hadronic contributions are substantially modified by low 
energy strong interaction effects, making the perturbatively obtained result unreliable. 
Therefore, we use two different improved schemes. In the first scheme, we start from the 
effective MS value at the Z pole, a = a eff (m z )|j^g ~ 1/127.7, where only the leptonic and 
light quark contributions are included. This results in the counter term |H] 



5e 
e 



(4 



7T 



e 

l 5 ~G 



O 2 

A + log4- 



/ ^ fJ f m=i 

+4 + log + ( A + log -^-) - 21 ( A + log 



(26) 



with Xf — m z V rrif < m z and x t = m t . N[, is the colour factor, N/j = 1,3 for (s)leptons 
and (s)quarks, respectively. A denotes the UV divergence factor, A = 2/e — 7 E + log47r. 
The second procedure we call the Gp scheme, with a = \/2s 2 v m 2 v G f / ^ — 1/132.5 and 



be 1 s w 
— = ~^Z AA - - — dZ 
e 2 2cw 



ZA 



1 . 

-Ar , 

2 ' 



(27) 



where Ar contains the full MSSM one-loop corrections [9] and the leading two-loop QCD 
corrections ^Uj- The leading contribution Act = I1^ A (0) — KeYl^ A (m 2 z )/m 2 z in Ar cancels 
the uncertainty in the transverse photon self-energy. Both schemes agree within a few per 
mill in the final results. 
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4.1.2 Higgs sector 



In the MSSM two complex Higgs doublets Hi and H2 are required. The ratio of the 
vacuum expectation values is defined by 



tan (3 = — . 

Vl 



(2* 



The renormalization of this mixing angle is by no means trivial [11 . In a definition by 
a specific physical process technical difficulties are introduced. A process independent 
renormalization on the other hand leads to gauge dependence and/or numerical instabili- 
ties. A convenient choice, is the condition, that the pseudo-scalar Higgs field A does not 
mix with the Z vector boson for on-shell momenta [T2*l I13j. 

5 tan f3 1 



Im 



Re U. A o z (m A o 



(29) 



tan (3 m z sin 2/3 

However, the result is gauge dependent and can lead to big corrections for large tan (3. It 
can be improved by taking only the divergent part of the A°Z self energy, which makes 
it numerically more stable and in the class of R% gauges £ independent at one-loop level. 
We call the first scale independent definition "on-shell" in contrast to the second "DR 
renormalization" at a certain scale Q. Both definitions and the translation between them 
are used in this work. 

4.1.3 SUSY sector 



For the renormalization of the chargino, neutralino and sfermion sector, we closely follow 
ref.|14|. We slightly extend the scheme for the general MSSM with complex phases. For 
other on-shell renormalization procedures, see e.g-JT^J. After rotation from the interaction 
to the mass eigenstate basis, we introduce again wave function and mass counter terms 



Xi 



(% + \5Z%P L + \bZf 3 P R )x 3 



m 



A 1 



m x, + Sm x> > 



f^(l + l5Zf)f. 



m 



f 



m 2 ~ + 5m 2 ~, hif = e L ,e R , 



(30) 
(31) 



where \ stands for both, charginos or neutralinos. These counter terms are related to the 
corresponding self-energies by 



5m 

sz, 



-Re [m* (^(mlJ + n£(m| 4 )) + H 



SL (m 2 



+ H 



SR( 2 
" L x. 



(32) 



Re 



SL, 



tSR, 



SL 'm 2 - 



+ nfKJ) 



Xi- 



Xi 



SL, 



"< ; 1I'/KJ + m^m^U^ml) + m^n^(m^) + m^% K (m\ 



Xj 



Xi- 



tSR, 
"Xj J - L ij 



Xj- 



(33) 



(34) 



10 



5Z£ = SZ^(L^R), 6Z§ = 6Z$(L R) , Cij = 2/{m%- m\) , (35) 

with 5Z = 5Z ± /5Z° and Uu(k 2 ) the chargino/neutralino self-energies, respectively. The 
corresponding counter terms for the sfermions are given by 

5m) = Re U f '(m 2 j) , 5Z f ~ = -Retl f (mf) , for f = e L , e R , z> e . (36) 

Furthermore we have to define the neutralino and chargino rotation matrices at one-loop 
level. A process and scale independent fixing seems appropriate. We define the counter 
terms in such a way that they cancel the rotation, induced by the antihermitian parts of 
the off-diagonal wave function corrections. 



5N i:i 



^£^-5Z^)N k0) (37) 

k=l 

1 4 1 4 

SU * = I E { 5 ^' R - ( S ^i R T) U *i > Wi = 4 E ( 5 ^' L - ( 5 ^' L )*) Vki ■ (38) 



fe=i fc=i 



In it is shown that this fixing eqs. (|37p and (|38|). calculated within the Feynman- 
t'Hooft gauge, can be regarded as a gauge independent one. The counter terms for the 
on-shell masses and the rotation matrices directly yield the counter terms for the mass 
matrix elements 



1 4 ~ 
SYn = ^ E N m N h Re 



2 

l,n=l 



4 

) + ™x? n SK ? ) + nff(m| ? ) + nff(m|o ; 



(39) 



1 4 ~ 
= 2 E Re 



(,n=l 



m^±nf,(m-±) + m-±n^(m 2 ±) + IEf (m-±) + n^(m?± x 



(40) 



Note that in the presence of complex mass parameters this renormalization automatically 
includes the one-loop definition of CP violating phases. 

For a proper one-loop calculation each involved Lagrangian parameter requires a clear 
definition, i.e. a unique counter term. Since not all entries in the neutralino and chargino 
mass matrices are independent parameters, in general the counter terms eqs. (|3*9*j) and 
(pfUj) cannot be interpreted as the counter terms to the Lagrangian parameters given in 
the tree-level mass matrices. For example the parameter M cannot be defined by the 
counter term 8Y22 to the neutralino mass matrix and in the same computation by the 
chargino mass matrix counter term 5Xn. The same holds holds for fi and the parameters 
already fixed in the gauge and Higgs sector. Therefore, we define M and \i to be the 
parameters in the chargino mass matrix M = X n , fi = X 22 and further M' = Yu 
throughout the calculation. For all other elements finite shifts have to be taken into 
account, e.g. Y 22 = M + SX n - 8Y 22 = M + AM. The UV finiteness of these shifts 
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is a nontrivial check of this method. The so obtained one-loop corrected mass matrices 
give after diagonalization the one-loop on-shell neutralino and chargino masses. The 
corresponding rotation matrices have the appropriate values for the counter terms eqs. (|37jl 
and (J38j) . This procedure of defining on-shell parameters can be simply extended to 
the sfermion sector. In our simplified case without left-right mixing, there are two free 
parameters for three masses. We fix Mr and Mg in such a way, that the selectron masses 
do not obtain one-loop corrections 

5M 2 Li = 5mj L + 5 (m 2 z cos 2/3 Q - sin 2 6 W ) \ , (41) 
5M| x = bm\ R + 5 (m 2 z cos 2/5 sin 2 9 W ) . (42) 
Due to a finite shift in , the sneutrino mass obtains the one-loop correction 

Ami = SM 2 Li + ^5(mlcos2(3) -8m\. (43) 



5 QED corrections 

The full one-loop corrections include diagrams with virtual photons attached to the tree- 
level diagrams. These contributions are IR divergent and regularized by an infinitesimal 
photon mass A. Concerning chargino and neutralino production these diagrams cannot be 
separated from the residual weak corrections in a gauge invariant and UV finite way. This 
can be traced back to the tree-level selectron and sneutrino t-channel diagrams, which 
introduce the charged current coupling g = e/s w . The same effect can be observed in 
the SM, e.g. for W pair production ^7j. The cross sections become IR finite and thus 
physically meaningful only by inclusion of real photon emission. For the calculation of 
the real photonic corrections, we use the so-called phase-space slicing method [TH|- The 
singular soft and collinear parts in the bremsstrahlung phase space are separated from 
the finite region. Both contributions can be written proportional to the tree-level cross- 
section, up to small terms of 0(AE / \/s) and 0(A6), and performed analytically. The 
collinear (soft) singularities are regularized by the electron (infinitesimal photon) mass 
and cancel the corresponding terms in the virtual corrections. In the soft photon area 
only photons up to a CM energy AE are considered 

da soft = _ da t rC e a J d^ L ftf__rt_ + Q K _ Q _k%\ 2 ^ 



4?r 2 Jk9<AE k® \Pikj hk-f k 2 k iy 

whereas = — 1/0 stands for the electric charge of the chargino/neutralino respectively. 
The result can be expressed in terms of the soft photon integrals given in ref. [7j 

da so ft = -da-c£_ ( 5 ISR + 5 fsr + S - FSR) ( (45) 
^soft = Ipwi + Ip2P2 ~~ 2-fpip 2 5 ^sof? = Q\ (-^fcifci + h 2 k 2 — 2Ik 1 k 2 ) ) (46) 

^soft = ^Qx (?Plki + Ip 2 k2 ~~ Ipik 2 ~ Ip 2 ki) ■ (47) 
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The collinear part takes hard photons in a small angle AO around the beam axis into 
account 



da coll ( Pl ,p 2 ,£-,£ + ) = f- 



a 



2tt 



l-2AE/y/s 







dx G a (s, x, AO) . 



a=± 

da tvcc (x Pl ,p 2 ,a^^ + ) + / da trcc ( Pl ,xp 2 ,t-,at 



with 



G+(s,x,A0) 



1 + x 2 

1 — X 



log 



( sAO 2 
\ Ami 



G-(s,x, AO) 



1 — x . 



(4£ 



(49) 



The finite hard bremsstrahlung has to be calculated by integration of the squared tree- 
level matrix-element for e + e~ — > XiXjl over the three-particle final-state phase-space. 
The complete 0(a) corrections can then be written as sum of virtual, soft, collinear, and 
finite contributions, all depending on unphysical auxiliary parameters. 

Aa 0{a ) = J (da virt (A) + da soft (A, AE)) + J da coll (A£, AO) + J da finitc (A£, AO) (50) 

Summing up the contributions in Aao{a), we obtain a cut-off independent result. This 
has been checked analytically for A and numerically for AE and AO in the intervals 
KT 5 < AE/^f7s < KT 2 and 10~ 3 < AO < 10~ 2 . 

For precise predictions of chargino and neutralino pair production higher orders beyond 
0(a) have to be taken into account. The structure function formalism ^| provides the 
possibility of defining process-independent logarithmic QED corrections, which originate 
from collinear virtual and real photons radiated off the incoming electron-positron beams. 



da 



tree 



da L 



dxJ dx 2 T^(x 1 ,Q 2 )T^(x 2 ,Q 2 ) liW'-i.ra, .r 2l , 
o Jo 



(51) 



We use the leading-log structure function up to 0(a 3 ), given in ref.[20j. 

exp(-§/3 7B + , 



x,Q 2 



r(i + f; 



8. , 8 2 ( . , , 21ogx 3, 

- 1(1 + x) + P - ( - 2(1 + x) log(l -x)- + -(1 



x) log X 



X 



+( 



1 / 9 7T 2 3 1 1 

"T^ 1 +x ) ( 32 ~ 12 + 4 log ^ ~ x "> + 2 log2 ^ ~ x )~ ^ l °Sx\og(l - x) 



H log 2 x Liofl 

16 & 4 v 



11 + X V 3 , 1 , 2 1, , H N \ 

2 l-x V 8 12 g x_ 2 gX g ^ V 



-^(1 -s)(log(l -x) + +^(5-3x)logx 



(52) 
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with the gamma function T, the Euler constant 7^ ~ 0.577216, and (3 = — (log% - 1). 
For simplicity we fix the free scale Q 2 = s throughout the calculations. Since 8 contains in 
addition to the log ^ a constant term also non-leading log terms are included to form the 
correct ISR soft-photon pole, see also eq. (|5T|. Furthermore, the soft-photon contributions 
are resummed up to all orders in perturbation theory. 

Since the above calculated correction Aaot a ) already contains the universal terms of 0(a) 



Jda™"' 1 = ^dx 1 ^dx 2 T^ 1 (x 1 ,Q 2 )T^ 1 (x 2 ,Q 2 ) Jda tiee (x lPe -,x 2 p e 



(53) 



with 



T^\x,Q 2 ) = §lim{5(l-x)[| + 21og( e )] + 0(l-x-6)i±^}, (54) 
4 e->o 2 1 — X 

we have to subtract them from the complete corrections to avoid double counting. 

J da complctc = J da tree + Aa 0{ a) + J (da univ - da"™' 1 ) (55) 



Definition of weak and QED corrections 

In spite of the impossibility to disentangle the different contributions to the 0(a) correc- 
tions, it is of special interest to distinguish the genuine weak corrections from the large 
and on experimental cuts dependent photon part. In fact, there is no unique way of doing 
this. A naive treatment would be to take the pure virtual corrections and set the photon 
mass equal to a typical scale of the corresponding process A = Q. However, this leaves 
us with enhanced Sudakov double-logarithms log 2 from virtual soft photons attached 
to the incoming beams, which are cancelled by the corresponding real, soft part. In the 
following we take the sum of virtual and soft corrections and extract the AE depen- 
dent terms as well as the contributions proportional to L e = log stemming from the 
collinear virtual+soft photons, eq. (|5*3|) 1 . 



da 



weak 



da 



virt+soft 



(AE) - -da 

7T 



tree 



AAE 2 
log (L e 



1 + A. 



2 



(56) 



The term A 7 takes the cut-off dependent terms from final state radiation (FSR) and 
ISR-FSR interference of eq. ()45|) into account. The sum da tree + da weak is identical to the 
"reduced genuine SUSY cross section" within the SPA convention The full corrected 
cross-section can now be obtained by the sum 



J d(T complete = f ^^tree + d(J wcak) + f d(J non-umv + J 



(57) 



1 The definitions in the previous work slightly differ from those given here. 
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with the non-universal QED corrections 



f d(J non-umv = ^(7 + (AS) — / da^' 1 + 

f d(7 tree 



O 

TT 



4Ai?^ 3 
log (L e - 1 + A 7 ) + -L e 

S — 



(58) 



A second way to highlight weak corrections is to compare the complete corrected cross 
section (including the hard photon radiation) with an improved tree-level d<j tree + ISR that 
already contains the universal QED corrections. 

J ^^.complete _ J ^^.trcc+ISR _|_ J ^^.residual (59) 
/ d(T tree+ISR = / d(T tr ° C + / d(T Univ , / d(T rcsidual = J ^weak + J d(J non-univ _ 

The advantage of this definition is that it does not require a more or less superficial split- 
ting of virtual and real corrections. On the other hand, the "residual" corrections include 
the non-universal QED corrections. They are in general small, but can be comparable to 
the loop corrections, especially the ISR-FSR terms. Furthermore, it can be inconvenient 
for technical reasons to include the hard bremsstrahlung process in the definition of a 
"weak" correction. 



6 Numerical results 

For the numerical analysis we concentrate especially on the SPSla' point, proposed in the 
SPA project jlj. It is close to the original Snowmass SPSla scenario and compatible with 
all available precision data and actual mass and cosmological bounds. The parameters in 
the SPA convention are defined in the DR scheme at the scale of Q = 1 TeV. A translation 
from these parameters to our on-shell definition can be simply performed by subtraction 
of the corresponding counter terms, i.e. V os = V(Q) — 5V(Q). The used values in this 
work can be seen in Table ^ For all other parameters that are free of renormalization 
conditions, the DR or on-shell values can be used. The difference is of higher order for 
the current processes. 

Since we use an on-shell renormalization, the appropriate tree-level for the one-loop cal- 
culation is given in terms of on-shell parameters. On the other hand, our original input 
are the DR parameters of the SPSla' scenario. We therefore show the corrections relative 
to the tree-level in the SPA convention, i.e. calculated in terms of on-shell masses and DR 
values for all couplings. Using this tree-level definition, the relative corrections are the 
same compared with those calculated in other renormalization schemes, up to terms of 
higher order. For this purpose the used DR value at Q=l TeV of the fine structure con- 
stant is a = 1/124.997. In the presented numerics, we use for the charge renormalization 
the a eS (m z )\jj^ scheme for neutralino production and the Gf scheme for chargino produc- 
tion, as discussed in section H.l .11 By calculating the same cross section in both schemes, 
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103.22 
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193.31 


197.03 
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tan (3 
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M El 


115.59 


117.71 


Ml. 


181.25 


183.98 



M 


OS 


M 


OS 


m xt 
m xf 

m e L 

m~ eR 


184.2 
421.1 
190.1 
125.2 
172.8 


m xl 


97.75 
184.4 
406.9 
419.5 



Table 1: Parameters of the SPSla' scenario in the DR and on-shell scheme and particle 
masses. 

we find good agreement within a few per-mill in the final results. The Figs. (J7J) and (JBJ) 
show the total cross sections for chargino and neutralino pair production at the tree-level 
in the SPA convention, and with weak and full corrections. Additionally, Fig. J7J) shows 
the complete corrections to the improved tree-level, where ISR is already included. Since 
the soft-photon pole is already absorbed into the tree-level, we have moderate corrections 
even near the threshold. 
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Figure 7: Left: Total cross sections for chargino pair production at tree-level {dotted}, 
with weak {dashed}, and with complete corrections {solid}. Right: Complete corrections 
relative to the improved tree-level above the particular production threshold. 



Comparing the various corrections to the total tree-level cross section for the different 
production channels Figs. (0 E3 ITT]) some common characteristics can be recognized. 
Near the threshold the negative soft photon contributions dominate. Far away from the 
threshold, the positive universal QED corrections partially cancel the large and negative 
weak contributions. The almost constant non-universal QED corrections are small and 
in comparison with other corrections often negligible. Due to the marginal non-universal 
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Figure 8: Total cross sections for neutralino pair production at tree-level {dotted}, with 
weak {dashed}, and with complete corrections {solid}. 




Figure 9: Radiative corrections for chargino pair production. The {full, dashed, dot- 
ted, dash-dotted} line corresponds to the {complete, weak, non-universal QED, universal 
QED} corrections to the total tree-level cross section. 



QED corrections, the differences between the two proposed ways to highlight "genuine 
weak" corrections are quite small. However, this does not have to be true any longer 
if cuts on the phase space are applied or distributions in particle energies or scattering 
angles are discussed. Moreover, we present results for the forward-backward and left- 
right asymmetry A-fb arid Alr of chargino production in the Figs. (JT21 fT3j) . These are 
defined by 

cr F -o- B o-l-ctr 

-4fb = ■ , A LR = ■ , (61) 

a F + cr B vl + cr 

where a-p = <x(cos > 0) with 6L ^ is the angle between the incoming electron and the 
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Figure 10: Radiative corrections for chargino and neutralino pair production. The {full, 
dashed, dotted, dash-dotted} line corresponds to the {complete, weak, non-universal 
QED, universal QED} corrections to the total tree-level cross section. 




Figure 11: Radiative corrections for neutralino pair production. The {full, dashed, dot- 
ted, dash-dotted} line corresponds to the {complete, weak, non-universal QED, universal 
QED} corrections to the total tree-level cross section. 



outgoing Xi in the CMS frame. cr L / R denotes the total cross section for left /right-handed 
electrons and unpolarized positrons. Both asymmetries obtain sizeable corrections for all 
three production channels. The kinks in the lines can be traced back to so-called normal 
and anomalous thresholds. 
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Figure 12: Left: Forward-backward asymmetry for chargino pair production at tree-level 
{dotted}, with weak {dashed}, and with complete corrections {solid}. Right: Complete 
corrections to A FB . 
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Figure 13: Complete corrections to Alr for chargino pair production. 



7 Conclusions 



We have presented in detail the calculation of 0(a) radiative corrections to the pair 
production of charginos and neutralinos within the Minimal Supersymmetric Standard 
Model. We discussed a possible separation of weak and QED corrections. Although 
only the sum of virtual and real bremsstrahlung corrections have physical meaning, a 
separation of these two contributions has its advantages. The QED corrections are treated 
numerically by integration over the phase space with an additional photon in the final 
state. Quite often a cancellation between QED and weak corrections occurs leading to a 
small total correction. Furthermore, we used the structure function formalism to include 



19 



higher order QED corrections and to split the QED corrections into a universal and a 
non-universal part. The non-universal QED corrections contain the integration over the 
three particle final state and are in the case of the total cross-section often negligible. 
The numerical analysis is based on the SPSla' benchmark scenario, proposed within the 
SPA project. The DR parameters of the SPA convention are translated into a set of 
on-shell parameters, which serve as input for our on-shell renormalization scheme. Such 
translation tools make it possible to compare cross sections, masses, etc. calculated in 
different renormalization schemes. The presented numerics shows results for the total 
cross-section, left-right, and forward-backward asymmetries. In all three cases we have 
sizeable radiative corrections around the 10% and in some cases even higher. 
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